Anderson localization of solitons 
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At low temperature, a quasi-one-dimensional ensemble of atoms with attractive interaction forms 
a bright soliton. When exposed to a weak and smooth external potential, the shape of the soliton 
is hardly modified, but its center-of-mass motion is affected. We show that in a spatially correlated 
disordered potential, the quantum motion of a bright soliton displays Anderson localization. The 
localization length can be much larger than the soliton size and could be observed experimentally. 
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At zero temperature, cold atoms interacting attrac- 
tively in a one dimensional (ID) system tend to cluster 
together, forming a bright soliton. Explicit solutions of 
the many-body problem can be found in some cases, for 
example for contact interactions Using external po- 
tentials, it has been experimentally shown how to put 
solitons in motion Q ■ What happens to a soliton exposed 
to a disordered potential? If the potential is strong, it will 
destroy the soliton. If it is sufficiently weak and smooth 
not to perturb the soliton shape, one expects the soli- 
ton to undergo multiple scattering, diffusive motion and 
possibly Anderson localization [3j. Indeed, propagation 
of waves in a disordered potential is profoundly affected 
by Anderson localization. Multiple scattering on random 
defects yields exponentially localized density profiles and 
a suppression of the usual diffusive transport associated 
with incoherent wave scattering 0]. In ID Anderson lo- 
calization is a ubiquitous phenomenon [B[ , which has been 
recently observed for cold atomic matter waves [6[ . It is 
important to understand how it is modified when inter- 
actions between particles are taken into account. 

We consider a Bose-Einstcin condensate in a quasi-ID 
geometry. Within mean-field theory, it is described by 
the Gross-Pitaevskii energy functional 
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in units of Eq = AmuP^a 2 , Iq = h/2\a\mu}±, and to = 
h/4a 2 mui' 2 _ for energy, length and time, respectively. 
Here, ui±_ denotes the transverse harmonic confinement 
frequency, a the atomic s-wave scattering length, and (i 
the chemical potential. The cases of repulsive and at- 
tractive atomic interaction are covered by g = ±1. 

The dynamics is extremely different in both cases, so 
that we first discuss the case of attractive interaction, 
g = — 1 . The ground state of ((T|) is the bright soliton 0] 
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normalized to the total number of particles N. The chem- 
ical potential is fi = —N 2 /8 and the soliton width is 
£ = 2/N. This ground-state solution has an arbitrary 
center-of-mass (CM) position q and an arbitrary global 
phase 9 that spontaneously break the translational and 
the U(l) gauge symmetry of the energy functional (JTJ) , 
respectively. These degrees of freedom appear as zero- 
energy modes of Bogoliubov theory, and their quantum 
dynamics requires special attention 0, 0] . 

The energy functional ([I]) is no longer translation in- 
variant when a potential term J dzV(z)\4>\ 2 is added. If 
V(z) is sufficiently weak and smooth, the soliton shape 
remains unchanged to lowest order in V, and only its CM 
position q is affected. In this Letter, we show that the 
quantum dynamics of q in a disorder potential leads to 
Anderson localization of the soliton, over a localization 
length possibly longer than the soliton size. 

An intuitive picture emerges using a simple ansatz in 
terms of collective coordinates (cf. [10j in the context of 
Bloch oscillations). Inserting <f>(z\ q, P q ) = e qPqZ / N 4>q(z — 
q) into ([T]) leads to an effective quantum Hamiltonian 
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It describes the entire soliton as an object of mass iV" 
evolving in an effective potential / dzV (z)\(f>o(q — z)\ 2 
that is the convolution of the bare potential with the 
soliton density. 

This simple ansatz yields no information on the re- 
maining degrees of freedom. Therefore, we apply a more 
complete analysis, expanding the energy functional ((T|) 
to second order in deviations from the ground-state so- 
lution ([2]), as in Q, e.g. Diagonalization of the quadratic 
Hamiltonian results in the eigenvalue problem for the 
non-hcrmitian operator 
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Its right eigenvectors (u n ,v n ) and corresponding adjoint 



modes (u: 
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) build a basis that spans the functional 
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space of ((f), </>*). For all non-zero eigenvalues E n , the 
adjoint modes are left eigenvectors of C. This is no 
longer true for the zero-energy modes. The first zero 
mode (u e ,vg) = idg((j) , is related to the global U(l) 
gauge invariance 4>e~ ie — ► 0e~ l( - 6 ' +e * 1 broken by the clas- 
sical solution © and its adjoint mode is well known, 
(uf,vf) = d N {(j) 0l (j)l) @, Q. The other zero mode, 
(u q ,Vq) = id q {4>Q, 4>q), originates from the translational 
invariance q — > q + e also broken by the soliton solu- 
tion. To find the adjoint mode, one solves £(w^ d ,i;^ d ) = 
M _1 (u q , v q ), where M is determined by the requirement 
(u q d \u q ) - (v q d \v q ) = 1. This ensures that (u^vt*^ 1 ) is 
orthogonal to all eigenvectors of C with E n ^ [a, |9( . It 
is easy to verify that 
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and M = N, the mass of the system. 

Following Dziarmaga Q we may now perform an ex- 
pansion around the ground state, 



4> = 0o + Penf + Pquf + 



E 



(6„u„ + b* n v* n ) , (6) 



where all modes are implicit functions of 9 and q. In- 
serting ([B]) in the energy functional |T]) and expanding to 
second order in Pg , P q and &„ (requiring P 9 £ <C AT for the 
CM momentum) results in the Bogoliubov Hamiltonian, 
whose quantum version reads 
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where Pg = N — N = —idg and P q = —id q . 

Because 9 can take non-perturbatively large values, we 
can work in a subspacc of Hilbcrt space with definite 
N = N. The minus sign in front of Pg in ((?]) arises 
because the bright soliton © is a saddle point of the 
energy functional (fTJ. It has no consequences since the 
number of particles is fixed. 

When the energy functional is supplemented with a po- 
tential term J dzV(z)\4>\ 2 , this perturbation commutes 
with N, so that we can focus on its incidence on the 
second and third term in (J7]), i.e., the CM motion and 
Bogoliubov excitations. A smooth external potential can 
only slightly distort the soliton. If the potential energy 
variation across the soliton is of the order of the variance 
| Vb| of the disordered potential, then \V \ < = N 2 /8 
is a sufficient condition for the soliton shape to be only 
weakly modified by the external potential [ll|. Moreover, 
the large energy gap \/j,\ above the ground state [l2| makes 
Bogoliubov excitations by the potential negligible. The 
soliton shape thus follows adiabatically and reversibly, 
i.e. without heating, the variations of the external po- 
tential. Then, the only degree of freedom affected by 
the external potential is the soliton position q. Inserting 



the expansion ([6]) into the energy functional, expanding 
up to quadratic terms, and quantizing, we arrive at the 
effective Hamiltonian ([3]) to leading order in l/N and V. 

In the following, we study the case of A^ = 100 Li 7 
atoms with scattering length a = — 3nm in a trans- 
verse harmonic trap with lu±_ = 2ix ■ 5 kHz. Then, units 
for energy, length, and time are Eq = 1.28 • 10~ i Huj±, 
Iq = 47.8 iim and to = 0.25 s. A soliton ([2|) of size 
£ = 2/7V = 0.02 (« 1 /im) is initially prepared in a large 
axial harmonic trap with uj z = 100 (27r • 64 Hz). The 
trapping potential is sufficiently small not to distort the 
soliton, whose CM occupies the ground state of a har- 
monic oscillator with frequency lj z . When the trap is 
turned off, the soliton position starts its quantum dy- 
namics with the corresponding momentum distribution 



7T (fc) 



exp [-k 2 /Ak 2 ] , Afc 2 = Nu z , (8) 



and begins to explore the disordered potential. 

In ID random potentials, Anderson localization is 
generic: the amplitude of every plane wave with wave 
vector k decreases asymptotically as exp{— j(k)\q\/ 2}. 
The inverse localization length j(k) can be calculated 
analytically in the weak-disorder limit (see below) [l3j |. 
Equivalently, the energy spectrum is discrete and dense 
(in the limit of infinitely large systems) with exponen- 
tially localized eigenstates. 

Disorder potentials are completely characterized by 
correlation functions V{z{) . . . V(z n ) where the overbar 
denotes an ensemble average over disorder realizations. 
The average potential value shifts the origin of energy 
and can always be set to zero, V(z) = 0. The pair cor- 
relator can be written as V(z')V(z' + z) = V 2 C(z/ao), 
where Vq measures the potential strength, and <jq the 
spatial correlation length. Higher-order correlations are 
required to fully describe non-Gaussian disorder such as 
the optical speckle potential considered in the following. 

Optical speckle yields a light-shift potential V(z) oc 
x|-E(z)| 2 proportional to the intensity of the light field 
E(z) and to the atomic polarizability x, whose sign de- 
pends on the detuning of the external light frequency 
from the atomic resonance. At fixed detuning, the po- 
tential features either random peaks (the "blue-detuned" 
case) or wells ("red-detuned"). The potential distribu- 
tion is asymmetric, and the importance of odd moments 
can be probed by comparing the blue- and red-detuned 
cases for fixed |V |. We use |V | = 8 • 10~ 5 H = 0.1 in 
the following. A ID speckle potential has the pair cor- 
relation function C(y) = [sin(y)/y] 2 , with a correlation 
length that can be as short as 0.26 ^m Q or cto = 0.0056 
in our units. In fc-space, the corresponding power spec- 
trum reads Vv{k) = ir<r V 2 (1 - \ka \/2) 6(1 - \ka \/2). 

Our Hamiltonian j3]) shows that the bright soliton sees 
a convoluted disorder potential whose fc-space compo- 
nents are the product (A%fc£/2)/ sinh(7rfc£/2) x Vk of 
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Figure 1: (Color online) Inverse localization length of a 
bright soliton versus its momentum k in units of its width 
£ = 0.02 (~ 1 jtim) in a speckle potential with correlation 
length (Jo = 0.28£. Circles are numerical results obtained by 
diagonalizing the effective Hamiltonian Q, solid lines are nu- 
merical results from a transfer-matrix calculation. The bare 
potential amplitude is Vb = +/ — 0.1 (blue/red curve) in our 
units, corresponding to ±8 ■ 10~ J |/x|. Inset: same data in log- 
arithmic scale. The lowest-order perturbation result, Eq. @, 
is shown as a black dashed line. 



given by eq. ([8]) . Each /c-component localizes at long time 
into the asymptotic form cxp{— 7(fc)|g|/2}, implying that 
the superposition will also localize. If the extension of the 
initial state is small compared to the typical localization 
length (which is the case here), one can neglect the phase 
correlations between the various fc-components at long 
time. The final density distribution, ensemble-averaged 
over disorder realizations, is then given by: 



dfc^exp{-7(fc)|g|}7ro(/c). 



W{q) = 



(12) 

Using (fTTj) this integral can be evaluated by a saddle- 
point argument from which it becomes apparent that at 
distance q, one finds the soliton with an initial momen- 
tum k q such that \q\ = l/ r y(k q ) provided aAk ^> 1. The 
probability distribution for the soliton position then takes 
the form 



W(q) 



y/2 exp[-fc2/Afc 2 
eevAfc \q\ 



(13) 



Fourier components from soliton density and speckle po- 
tential. To second order (Born approximation) in the po- 
tential strength, the inverse localization length [l3| then 
reads 



7 (fc) 



N±e^oV 2 q~\kvo\) 
[sinh(7rfc£)] 2 



6(1 - |fcro|). (9) 



For a short correlation range Co -C £, the fc-dependence 
due to the soliton convolution dominates and the bare 
speckle can be approximated by its white-noise limit 
7-V(0) = ttuoVq. The soliton width £ takes over as the 
new effective correlation length scale. For 1/(tt£) < k < 
1 /(To one can use the approximate expression 

7 (fc) « (2ttN 2 2 P v (0) e- 2 ^ k . (10) 

The lowest-order perturbation result © can be com- 
pared with exact data computed numerically, both by ex- 
act diagonalization and transfer matrix methods [171 ] for 
the Hamiltonian Q. Figure [1] confirms that Anderson lo- 
calization is observed for all k- values, but also shows that 
the localization lengths for the blue- and red-detuned po- 
tential differ by up to an order of magnitude for fc£ > 0.5. 
Consequently, they also differ from the lowest order per- 
turbation results. Thus, perturbation theory cannot be 
expected to apply, even though the disorder potential is 
much smaller than the kinetic energy of the soliton. But 
the salient feature here is the rapid exponential decrease 
of j(k), clearly visible in the inset of Figurc[T] Thus, the 
inverse localization lengths still obeys 



7 = 70 exp(— ak) 



(11) 



with a cx £, the natural length scale of the problem. 

Let us return to the CM dynamics. The initial state is 
the superposition of various fc-components, with weights 



where k q = ln( 7 o|g|)/(cv) from (jlip . In the regime of 
interest, this is an almost perfect algebraic decrease as 
|(7| _1 with a small logarithmic correction (22l |. 

These predictions have been tested by numerical in- 
tegration of the Schrodinger equation with the Hamilto- 
nian ([3]) starting from the initial Gaussian wavefunction 
corresponding to (J5J), and averaging over 250 disorder 
realizations. As shown in Fig. [3J the final probability 
distribution for the soliton position follows quite well the 
predicted algebraic decay. The time scale required to 
observe a stationary localized state around position q is 
given by r(q) = \q\/v(k q ) = N\q\/k q = a|g|/ln(7 |g<|). 

Finally, we turn to the case of repulsive interac- 
tions [HI by taking g = +1 in {!]). The BEC ground 
state in a harmonic trap extends over the Thomas-Fermi 
radius. How such a condensate wave-packet expands in a 
disordered potential has been studied for different dimen- 
sionalities [14, 
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161 ] - Because of repulsive interaction, 
the entire condensate is an extended object that requires 
a field-like description, in sharp contrast with a bright 
soliton that features particle-like properties. 

There exists, however, a dark soliton solution 4>q oc 
tanh[(a; — q)/£] where £ = l/y/p and p is the background 
density. Then, one can use the non-perturbative descrip- 
tion presented in 0] or a collective-coordinate method 
[l9| to obtain an effective Hamiltonian for the dynam- 
ics of the dark soliton in the presence of a weak disorder 
potential: 



p 2 
g 

2IMI 



\M\ 
4£ 



dz- 



V(z) 



[cosh((z - q)/t)Y 



(14) 



The mass of the dark soliton is negative, M = — 4p£, and 
its modulus equals twice the number of particles missing 
in the soliton notch. This effective Hamiltonian is valid 
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Figure 2: (Color online) Probability distribution of soliton 
position q in units of its width £ = 0.02 (w 1 ^tm). Blue and 
red lines: Data from numerical integration of the Schrodinger 
equation with the Hamiltonian © starting in the ground state 
of an axial harmonic trap with uj, — 100 (2n ■ 64 Hz) and 
lasting At = 20 (5 s) for a blue/red-detuned speckle potential 
of amplitude Vo = ±8 • 10"V (2tt ■ 0.064 Hz). Data averaged 
over 250 different realizations of the disorder potential. Black 
dashed line: saddle point result (|13p for parameters in the 
Born approximation. The full integral (|12[) yields the same 
curve for \q\ S> £. The double-logarithmic inset plot shows 
the dominantly algebraic decrease as \q\~ . 



only when the velocity of the CM of the soliton is much 
smaller than the sound velocity c = ^fp in the conden- 
sate. At velocity comparable to c, the shape of the soliton 
changes, making the analysis more difficult. 

Equation (fl4|) shows that we obtain the same form of 
the effective Hamiltonian as for the bright soliton. At 
first sight, it predicts similar Anderson localization ef- 
fects. However, there is a fundamental difference between 
dark and bright solitons. The bright soliton ((2]), being 
the ground state of the N particle system, is protected 
by a large gap from quasi-particle excitations. In that 
respect, entire solitons can undergo quantum dynamics 
and localize like Rubidium atoms [6j . For the dark soliton 
a gap can only be imposed by the boundary conditions 
in a finite-size system, and is inversely proportional to 
the system size, making radiation of Bogoliubov excita- 
tions [2o| and decoherence much easier. Furthermore, the 
dark soliton is an excited state, as signaled by the global 
minus sign of the effective Hamiltonian p4[) . Interactions 
with thermal cloud may accelerate the soliton and make 
it disappear 2l| . Estimating the relevant time scale and 
evaluating its effect on the localization dynamics is left 
for future work. 

In summary we have shown that the center of mass of 
a bright soliton may undergo Anderson localization in a 
smooth disorder potential. In realistic situations where 
the soliton wavepacket is prepared in a small region of 
space, this leads to an essentially algebraic localization 
of the ensemble-averaged atomic density at long times. 
We emphasize that the effects discussed here are be- 
yond standard mean-field description: while the soliton's 
shape is described by a mean-field theory, its center of 



mass is treated quantum mechanically. This has impor- 
tant consequences: although the one-body density matrix 
will display Anderson localization as shown in Fig. [2j a 
single realization of the experiment is expected to find a 
single soliton at a given random position, with a proba- 
bility density given by W(q), eq. (TT3"]) . 
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